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FIELD OF THE INVENTION 

The invention relates to non-invasive in-vivo methods and apparatus for determining 
the concentration of a substance in a body. 

5 BACKGROUND OF THE INVENTION 

Non-invasive methods for assaying a "target" analyte, such as for example glucose, 
comprised in a region of body tissue are known in tiie art. In a near infrared spectroscopy 
(NmS) method, light at a plurality of different wavelengtiis in a near infrared band of 
wavelengtiis is transmitted into a tissue region of tiie body to assay a target analyte in the 

10 tissue region. Ught at at least one of tiie wavelengtiis, a "target wavelengfli" is absorbed or 
scattered by tiie target analyte. Intensity of light at tiie different wavelengtiis tiiat is transmitted 
through tfie tissue region or scattered out of tiie tissue region is measured. The measured 
uitensities are used to isolate and determine tiie contribution of the target analyte to an 
absorption or scattering coefficient of flie tissue region at the target wavelengtii in tiie 

15 presence of contributions to the absorption or scattering coefficient by otiier "interfering" 
analytes m the tissue. Known values for the absorption or scattering cross sections of tiie 
target analyte at tiie target wavelengtii and tiie determined contiibution of tiie target analyte to 
the absorption or scattering component are used to assay the target component in tiie tissue. 

However, MRS methods provide concentration measurements of a target analyte in 

20 tissue tiiat are averages over relatively long optical patii lengths tiurough tiie tissue of light 
used to acquire the measurements. As a result, NIRS metiiods and technologies generally 
suffer from poor spatial resolution. In addition, NIRS signals tend to suffa- from noise 
generated by scattering of light at tissue interfaces, such as tiie skin, and tissue 
inhomogeneities. NIRS metiiods tend tiierefore to exhibit relatively poor signal to noise ratios. 

25 An article by G. Yoon, et al. "Determination Of Glucose Concenti^tion in a Scattering 

Medium Based on Selected Wavelengtiis by Use Of an Overtone Absorption Band", in 
APPLIED OPTICS 1 March 2002; Vol. 41, No 7 describes an NIRS metiiod and device for 
assaying glucose in a tissue medium. The metiiod describes criteria for choosing discrete 
wavelengtiis for light used in performing an NIRS assay so as to reduce mfluence of 

30 interfering analytes in the tissue on flie results of tiie assay. Whereas tiie metiiod is based on 
measuring NIRS absorption spectira of tiie tissue medium, tiie tissue medium is assumed to be 
scattering as well as absorbing. The article describes a device for assaying glucose having a 
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light source and a detector that are used to measure absorption spectra for the medium that 
have their relative positions optimized so that "measured spectra can be independent of 
medium scattering". 

For many medical procedures it is advantageous to accurately determine concentration 
5 of a target analyte for tissue regions that are relatively spatiaUy localized. For example, in 
assaying glucose levels for a patient it is generally advantageous to measure glucose levels in 
blood. To acquire such measurements, the measurements should be spatially localized to a 
blood vessel or blood vessels so that the measurements are not "diluted", for example, by 
glucose levels in interstitial fluids. NIRS methods and devices, because of their relatively poor 

1 0 spatial resolution generally cannot provide such localized assays. 

Methods of measuring concentration of a target analyte in a tissue region using a time- 
resolved photoacoustic effect or optical coherence tomography (OCT) can provide 
measurements resolved to a relatively high spatial resolution. 

In a method using a time resolved photoacoustic effect, light at at least one wavelength 

1 5 for which light is absorbed or scattered by the target analyte is used to generate photoacoustic 
waves in the tissue region. Pressure produced by acoustic energy from the- photoacoustic 
waves that arrives at a suitable acoustic transducer or transducers is used to assay the target 
analyte at locations in the region at which the photoacoustic waves are generated. Locations at 
which the photoacoustic waves are generated can be determined to within about 10 microns 

20 axially, along a direction of propagation of the waves, and to within about 200 microns 
laterally. As a result, an assay of the target analyte can be spatially localized to relatively small 
volumes having axial dimensions of about 10 microns and lateral dimensions of about 200 
microns. 

In OCT, light from a semi-coherent light source comprised in an interferometer is split 
25 into a reference light beam and a light beam that aiuminates the tissue region. Light from tiie 
reference beam is reflected from a mirror to an "interference region" in the interferometer 
where it interferes with light scattered from the tissue region that reaches the interference 
region. An interference signal in the interference region is generated substantially only for 
reference and scattered light that reach the interference region after traveling substantially 
30 equal optical path lengths. As a result, flie interference signal is generated substantially only 
for light scattered from material located in a small volume of the tissue region for which the 
optical path lengths of scattered light and reference light are substantially equal. The 
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amplitude of the interference signal is substantially proportional to a scattering coefficient for 
material in the small volume and is used to assay the target analyte in the small volume. 
Optical coherence tomography can provide axial spatial resolution of about a micron and 
lateral resolution of about 3 microns. Spatial resolution of an assay provided by OCT assaying 
is therefore on the order of a small number of microns. 

However, absorption and scattering cross sections of a target and interfering analytes 
in a tissue region contribute to photoacoustic or OCT signals used to assay the target analyte. 
Accuracy of tiie assay is generally conq)romised if contributions to the signals from scattering 
cross sections of the analytes are not assessed and distinguished from contributions to die 
signals from absorption cross sections of the analytes. Prior art has not provided methods for 
assaying a target analyte in a tissue region responsive to photoacoustic or OCT signals for 
which scattering cross section contributions to the signals are assessed and distinguished from 
absorption cross section contributions to the signals. 

SUMMARY OF THE INVENTION 

An aspect of some embodiments of the present invention relates to providing assay 
apparatus that uses flie photoacoustic effect to assay a target analyte in a spatially localized 
tissue region and accounts for scattering of light used to generate the photoacoustic effect in 
determining the assay. 

An aspect of some embodiments of the present invention relates to providing assay 
apparatus that uses OCT signals to assay a target analyte in a spatially localized tissue region 
and accounts for scattering of light used to generate the OCT signals in determinmg the assay. 

An aspect of some embodiments of the present invention relates to providing a method 
for incorporating tiie effects of scattering of light in assaying a target analyte in a tissue region 
using the photoacoustic effect and/or OCT. 

An assay apparatus in accordance with an embodiment of the present invention 
comprises at least one light source that iUuminates the tissue region with light at each of a 
pluraUty of different wavelengths, hereinafter referred to as "mensuration wavelengths". For at 
least one of the mensuration wavelengths light is absorbed and/or scattered, optionally 
strongly, by the target analyte. 

In some embodiments of the invention, the assay apparatus comprises at least one 
acoustic transducer, which senses pressure in photoacoustic waves generated at different 
locations in tiie tissue region responsive to the light Alternatively or additionally, tiie assay 
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apparatus comprises an OCT interferometer. Interference signals are generated by the 
interferometer between light scattered from material at different locations in the tissue region 
and a reference beam of light provided by the ligjit source. An extinction coefficient for each 
mensuration wavelength is determined for the tissue region either responsive to signals 
5 generated by the at least one acoustic transducer or interference signals generated by the 
interferometer. 

The mensuration wavelengths are determined so that each extinction coefficient is 
dependent on concentration of at least one of a same plurality of "mensuration" analytes. One 
of the mensuration analytes is the target analyte and the remaining mensuration analytes are 

10 "interfering" analytes. Each extinction coefficient therefore defines an equation having as an 
unknown variable a concentration of at least one of the mensuration analytes. Together, the 
extinction coefficients define a plurality of simxjitaneous equations having as unknown 
variables concentrations of the mensuration analytes in the tissue region. The number of the 
plurality of mensuration wavelengths and therefore the number of simultaneous equations is 

15 equal to or greater than the number of the plurality of mensuration analytes. 

In accordance with an embodiment of the present invention, a wavelength dependent 
function, hereinafter a "scattering coefficient function", which is parameterized by at least one 
characteristic parameter, is used to provide a value for the scattering coefficient for at least 
one of the mensuration wavelengths fliat contributes to the extinction coefficient at the 

20 wavelength. The equation defined by the extinction coefficient for a given mensiiration 
wavelength comprises a term which is the scattering coefficient function evaluated at the 
mensuration wavelength. An assay of the target analyte is provided responsive to a solution of 
the simultaneous equations. 

In some embodiments of tiie present invention, at least one characteristic parameter of 
25 the scattering coefficient function is determined from information extraneous to information 
used to determine the set of simultaneous equations. In some embodiments of the present 
invention, at least one characteristic parameter of the scattering coefficient function is 
determined from an extinction coefficient determined from signals provided by the at least 
one acoustic transducer or altematively by the interferometer. 
30 In some embodiments of the present invention, the number of the plurality of 

mensuration wavelengths and therefore extinction coefficients is greater than the number of 
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the plurality of mensuration analytes and the simultaneous equations are used to determine at 
least one characteristic parameter of the scattering function. 

In some embodunraits of the invention, the scattering coefficient function is 
determined assuming that optical scattering in the tissue region is Mie scattering. 
5 hi some embodiments of tiie invention, the target analyte is glucose and an assay 

apparatus is used to provide in vivo measurements of glucose in a blood vessel in the body of 
a patient. 

There is therefore provided in accordance with an embodiment of the present 
invention apparatus for assaying a target analyte in a localized tissue region fliat may include 

10 tile target and otiier analytes comprising: a light source tiiat illuminates tiie region with light at 
each of a pluraUty of wavelengths at which light is absorbed and/or scattered by tissue in tiie 
region wherein light at at least one of tiie wavelengths is absorbed or scattered by flie target 
analyte; a signal generator that generates signals responsive to intensity of the light from the 
light source at different locations in tiie localized region; and a controller tiiat: receives tiie 

15 generated signals; processes tiie signals to determine an extinction coefficient for light in flie 
localized region at each wavelengtii; and determines tiie concentiation of the target analyte 
responsive to a solution of a set of simultaneous equations having as unknown variables 
concentrations of a pluraUly of analytes in the region, one of which is tiie target analyte, 
wherein each equation m tiie set expresses a relationship between tiie extinction coefficient, 

20 tiie absorption coefficient and/or tiie reduced scattering coefficient for light at a different one 
of tiie plurality of wavelengths and at least one of tiie equations expresses a relationship 
between flie extinction coefficient and the reduced scattering coefficient. 

Optionally, tiie at least one equation that expresses a relationship between flie 
extinction coefficient and tiie reduced scattering coefficient includes a dependence on the 

25 absorption coefficient. 

AdditionaUy or alternatively tiie reduced scattering coefficient at at least one of flie 
wavelengths is a measured value of the reduced scattering coefficient. 

hi some embodiments of the present invention, the reduced scattering coefficient at at 
least one of the wavelengflis is a value determined responsive to an analytic expression. 
30 hi some embodiments of tiie present invention, tiie reduced scattering coefficient at at 

least one of flie wavelengths is expressed as an analytic function. Optionally, the analytic 
expression is a function of at least one unknown variable having a value determinable 
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responsive to a solution of the simultaneous equations. Optionally, the at least one unknown 
variable is a concentration of at least one of the target analyte and the other analytes. 

hi some embodiments of the present invention, the function comprises an expression 
of the form BX'^ where X represents the wavelength and B and C are constants. 
5 In some embodiments of the present invention, the signal generator comprises at least 

one acoustic transducer that generates signals responsive to acoustic energy that reaches the 
transducer from photoacoustic waves generated in the region by the light. 

In some embodiments of the present invention, the signal generator comprises an 
optical coherence tomography device flxat receives light from the light source that is scattered 
10 from the region and generates an interference signal responsive to an interference pattern 
between the scattered light and light from the light source reflected by a reflector. 

In some embodiments of the present invention, the controller identifies and locates the 
localized region in a larger region comprising the localized region. 

Optionally, to identify and locate the localized region the controller: controls the light 
15 source to illuminate the larger region with light that is absorbed by a component characteristic 
of the localized region; receives signals generated by the signal generator responsive to 
intensity of the light from the light source in different locations in the larger region; uses the 
signals to assay the characteristic component in different localized regions in the larger region; 
and identifies and locates the localized region responsive to the assay. 
20 Optionally, the apparatus comprises at least one acoustic transducer controllable to 

transmit ultrasound, and to identify and locate the localized region the controller: controls the 
at least one transducer to transmit ultrasound into the larger region; receives signals generated 
by the at least one acoustic transducer responsive to acoustic energy reflected by features in 
die larger region from the transmitted ultrasound; and uses the signals to identify and locate 
25 the features and thereby the localized region. 

In some embodiments of the present invention, the localized region is a bolus of blood. 
There is further provided in accordance witfi the present invention a method of 
assaying a target analyte in a region of body tissue that may include the target and other 
analytes comprising: determining an extinction coefficient for light at each of a plurality of 
30 different wavelengths at which light is absorbed and/or scattered by tissue in the region and 
wherein light at at least one of the wavelengths is absorbed and/or scattered by the analyte; 
providing a value or an analytic expression for the reduced scattering coefficient at each 
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wavelength; and determining the concentration of the target analyte responsive to a solution of 
a set of simultaneous equations having as unknown variables concentrations of a plurality of 
analytes in the region, one of which is the target analyte, wherein each equation in the set 
expresses a relationship between the extinction coefficient, the absorption coefficient and/or 
5 the reduced scattering coefficient for light at a different one of the plurality of wavelengths 
and at least one of the equations expresses a relationship between the extinction coefficient 
and the reduced scattering coefficient. 

Optionally, the at least one equation that expresses a relationship between the 
extinction coefficient and flie reduced scattering coefficient includes a dependence on the 
1 0 absorption coefficient 

Additionally or alternatively determining the extinction coefficient at at least one of 
the wavelengths of the plurality of wavelengths optionally comprises: fi-om a given location 
iUummating the region with light at the wavelength so as to generate photoacoustic waves in 
the region; determining a rate of decrease amplitude of llie generated photoacoustic waves 

15 with increase of distance in the tissue region firom the given location; and determining the 
extinction coefficient firom the determined rate of decrease. 

In some embodiments of the present invention, determining the extinction coefficient 
at at least one of tiie wavelengths of flie plurality of wavelengtiis comprises: from a given 
location illuminatmg the region with light at the wavelength; using optical coherence 

20 tomography to determine a rate of decrease of intensity of the light with increase of distance in 
the tissue region from the given location; and determining the extinction coefficient from the 
determined rate of decrease. 

In some embodiments of tiie present invention, the reduced scattering coefficient at at 
least one of the wavelengths is a measured value of the reduced scattering coefficient 
25 In some embodunents of the present invention, the reduced scattering coefficient at 

least one of the wavelengths is a value determined responsive to an analytic expression. 

In some embodiments of the present invention, the method comprises expressing the 
reduced scattering coefficient in at least one of the equations as an analytic fimction. 

In some embodiments of the present invention, die analytic expression is a function of 
30 at least one unknown variable having a value determinable responsive to a solution of the 
simultaneous equations. Optionally, tiie at least one unknown variable is a concentration of at 
least one of the target analyte and other analytes. 
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In some embodiments of the present invention, the analytic expression comprises an 
expression of the form BX-C where X represents the wavelengfli and B and C are constants. 

In some embodiments of the present invention, the method comprises identifying and 
locating the localized region in a larger region comprising the localized region. 

Optionally, identifying and locating the localized region comprises: illuminating the 
larger region with light that is absorbed by a component characteristic of the localized region; 
generating signals responsive to intensity of the light at different locations in the larger region; 
using the signals to assay the characteristic component in different localized regions in the 
larger region; and identifying and locating tilie localized region responsive to the assay. 

Additionally or alternatively, identifying and locating the localized region comprises: 
transmitting ultrasound into the larger region; generating signals responsive to acoustic energy 
reflected by features in the larger region from the transmitted ultrasound; and using the signals 
to identify and locate the features; using the identities and locations of the features to identify 
and locate the localized region. 

In some embodiments of the present invention, the localized region is a bolus of blood. 

BRIEF DESCRIPTION OF FIGURES 

Non-limiting examples of embodiments of the present invention are described below 
with ref»ence to figures attached hereto, which are listed following this paragraph. In the 
figures, identical structures, elements or parts that appear in more than one figure are generally 
labeled with a same numeral in all the figures in which they appear. Dimensions of 
components and features shown in the figures are chosen for convenience and clarity of 
presentation and are not necessarily shown to scale. 

Fig. 1 schematically shows an assay apparatus assaying glucose using the 
photoacoustic effect, in accordance with an embodiment of the present invention; and 

Fig. 2 schematically shows an assay apparatus that uses both the photoacoustic effect 
and OCT to assay glucose, in accordance with an embodiment of the present invention. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Fig. 1 schematically shows an assay apparatus 20, hereinafter referred to as a 
"glucometer", assaying glucose in a "target region" 22 of a body part 24 of a patient, in 
accordance with an embodiment of tiie invention. Target region 22 is optionaUy located in a 
region 26 of soft tissue of body part 24 and comprises a body fluid, such as for example 
interstitial fluid, having a concentiiation of glucose. Optionally, target region 22 is a volume of 
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body fluid having a concentration of glucose and region 26 is a region of a fluid cavity 
containing the body fluid. For example, as in Fig. 1, target region 22 is a bolus of blood and 
the fluid cavity a blood vessel 23. 

Glucometer 20 optionally comprises a controller 32, a light source 34, optionally 
5 located in the controller, and an optic fiber 36 coupled to the light source. An end 38 of fiber 
36 is optionally mounted to a support structure 40, hereinafter a "probe head", to which an 
acoustic transducer or array of transducers is mounted. Any of various appropriate acoustic 
transducers or array of transducers may be used in the practice of the invention. By way of 
example, in Fig. 1 probe head 40 has an array of acoustic transducers 42 positioned 
10 circumferentially around end 38 of optic fiber 36. Only two transducers of the array are 
shown. Probe head 40 is pressed to skin 44 of body part 24 to position end 38 of fiber 36 
close to or contiguous witii the body part and to acoustically couple acoustic transducers 42 to 
the body part. 

To assay glucose in blood bolus 22 conti^oller 32, optionaUy first, controls glucometer 
15 20 to locate blood vessel 23 and the bolus using any of varioiis metiiods known in the art, such 
as metiiods described in PCX publication WO 02/15776, the disclosure of which is 
incorporated herem by reference. For example, to locate blood bolus 22 controller 32 may 
control transducers 42 to radiate ultrasound into region 26. Contioller 32 processes signals 
generated by tiansducers 42 responsive to reflections of tiie radiated ultrasound from 
20 stiiictures in region 26 to determine location of blood vessel 23. Alternatively, controller 32 
may control light source 34 to illuminate tissue region 26 with light fliat is relatively strongly 
absorbed by blood. Since flie light is strongly absorbed by blood, photoacoustic waves are 
preferentially generated in blood vessel 23. Controller 32 processes signals generated by 
transducers 42 responsive to acoustic energy firom the photoacoustic waves to image features 
25 in region 26 and locate blood vessel 23 . 

In accordance with an embodunent of the invention, controller 32 then controls light 
source 34 to illuminate region 26 with at least one light pulse, represented by wavy arrows 50, 
at each of a plurality of mensuration wavelengtiis ^. The index i indicates a particular one 
of tiie Nx mensuration wavelengtiis and satisfies the condition 1 ^ i< N^,. 
30 The at least one pulse of light 50 at target wavelengfli "ki stimulates photoacoustic 

waves, schematically represented by starbursts 52, in tissue region 26 and bolus 22. 
Transducers 42 generate signals responsive to pressure in acoustic energy from photoacoustic 
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waves 52 that reach the transducers. The signals are transmitted to controUer 32. which 
processes the signals m accordance with an embodiment of the invention, as described below, 
to determine glucose concentration in target region 22. 

In some embodiments of the invention, the at least one pulse of light 50 transmitted at 
different mensuration wavelengths ^ is transmitted at different tunes to iUummate bolus 22. 
Li some embodiments of the invention the at least one pulse 50 comprises a train of pulses. In 
some embodiments of the invention the pulses in the train of Ught pulses at different 
mensuration wavelengths Xi are transmitted at different pulse repetition rates. Optionally, light 
pulse trains at different mensuration wavelengths are transmitted simultaneously. Signals 
generated by acoustic transducers 42 responsive to photoacoustic waves 52 that are stimulated 
by light pulse trains at different mensuration wavelengths are distinguished using signal 
processing techniques known in the art, such as appropriate heterodyning and phase locking 
techniques. 

Let mtensity of a light pulse 50 transmitted at a mensuration wavelength 7^ into body 
part 24, at a distance d from end 38 of fiber 36 be represented by I(Xi.d). Assuming that d is 
larger than the mean free path for photons at wavelength Xi, I(Ai,d) may be written, 

(1) I(;ii,d) = Io(Ai)exp(-ae(Ai)d) 

where 0^0^ is an extinction coefficient in the tissue of the body part for light at wavelength 
Xi and Io(Ai) is intensity of light in the Ugjit pulse at end 38 of fiber 36. The extinction 
coefficient is a function of an absorption coefficient ota(Ai) and a scattering coefficient OgCAi) 
in the tissue for light at the wavelength 7^. Under the assumption of the diffiision 
approximation, the extinction coefficient may be written 

(2) Ch(M) = [3aa(A^(oca(^ + a'sCXj))] 1/2 
where, 

(3) a'siXO = (l-g)cxs(^. 

a's(Xj) is referred to as a reduced scattering coefficient and g is an anisotropy factor. 

Equation (2) may be rearranged to provide an expression for the absorption coefficient 

(4) OaCXO = 1/2 {-a's(A^[ a's(Xi)2+(4/3)ote(Xi)2] 1/2 } . 

At each wavelengdi Xi, flie absorption coefficient a^(Xi) may be expressed as a sum of 
absorption coefficients of analytes in region 24 that absorb light at the wavelength Xi. The 
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absorption coefficient of a given analyte is a product of an absorption cross-section of the 
analyte for Ught at wavelength 7^ and concentration of the analyte in the body. Let the 
absorption cross-section of a "j-th" analyte at wavelength Xi be represented by ajCXj) and its 
concentration in blood bolus 22 by xj. In accordance with an embodiment of the present 
invention the mensuration wavelengths are chosen so that at each mensuration wavelength 
Xi, substantially only at least one of a same plurality of "Na" mensuration analytes contributes 
to otaCAi). One of the Na mensuration analytes aj(Xi) is glucose, the target analyte, and one of 
the mensuration wavelengths Ai is a target wavelength corresponding to the target analyte 
glucose for which target wavelength Ught is, optionally, strongly absorbed by glucose. The 
absorption cross section aj(^ for j=l and the wavelength Xi for i=l are arbitrarily assigned to 
represent respectively the absorption cross section for the target analyte glucose and the 
corresponding target wavelength. The absorption coefficient at wavelength 7^ may therefore 
be written, 

Na 

(5) aa(Xi)= J^<yMx.= l/2{a's(Xi)+[a's(>a)2+(4/3)0Ce(Ai)2]l/2}. 

j ^ 

The mensuration wavelengths provide a set of linear equations of the form of 

equation (5) in the Na unknown concentrations xj (1 ^ j ^ Na). The equations can be solved 

for any and aU the mensuration analyte concentrations xj, and in particular for concentration 

XI of glucose in blood bolus 22 if N;^, > Na and for each mensuration wavelength Ai of the 

Nx wavelengths the extinction coefficient ccgCAi) and reduced scattering coefficient a'sC^i) are 
known. 

To determine the extinction coefficient OeCA^) for blood bolus 22 for each mensuration 
wavelengths Xi in accordance with an embodiment of the invention, controUer 32 processes 
signals that transducers 42 generate responsive to pressure produced by photoacoustic waves 
52 at tiie transducers. 

Pressure sensed by acoustic sensors 42 responsive to photoacoustic waves 52 is time 
dependent. Pressure sensed at a time "t" following a time at which a light pulse 50 illuminates 
body part 24 arises from photoacoustic waves generated at locations in the body part for 
which distance "d" from acoustic sensors 42 is substantially equal to vt, where v is the speed 
of sound. (The transmission time of the light is negligible.) Let the pressure sensed by acoustic 
sensors 42 at time t responsive to a pulse of light 50 at wavelength Xi that illuminates tissue 
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region 26 be represented by P(^,t). Then for photoacoustic waves generated at locations in 
tissue region 26 at a distance d from transducers 42, P 0<^t) can be written: 

(6) F(Xi,t) = P(Xi,d/v) = Kota(Xi)I(Xi,d). 
where K is a constant of proportionality. 
Using equation 1 , equation (6) may be rewritten, 

(7) P(Xi,t) = P(Ai,d/v) = Kcea(Ai) {Io(Ad)exp(-oee(Xi)d)} . 

From the time dependence of P(^,t), controller 32 determines which of the signals 
generated by transducers 42 responsive to P(Ai.t) are generated responsive to photoacoustic 
waves originating at distances d from fiber end 38 corresponding to locations in bolus 22. 
(Distances d that correspond to bolus 22 are known from the location of blood vessel 23, 
which was determined as noted above.) From the signals responsive to photoacoustic waves 
52 originating inside bolus 22 controller 32 determines values for F(Ki,d) for a pluraUty of 
locations in blood bolus 22 at different distances d from end 38. The controUer uses the 
determined values for P(^,d) and equation (8) to determine a value for ae(A^. OptionaUy the 
determined value for cXgCA^i) is a best fit value that optimizes the fit of equation (7) to the 
determined values for P{7^,d). 

To determine the reduced scattering coefficient a's(^ for mensuration wavelength 'ki, 
in some embodiments of the present invention, scattering of light in blood bolus 22 is 
measured at the wavelength. In some embodiments of the invention the scattering coefficient 
is determined from a wavelength dependent analytic fimction, i.e. a scattering coefficient 
function parameterized by at least one characteristic parameter. Optionally, the scattering 
coefficient function is determined assuming that scattering of light is substantially Mie 
scattering. As a result, as is known in the art, dependaice of a's(Ai) on wavelength may be 
approximated by an expression of the form, 

(8) a's(Ai) =BAi-C. 

In some embodiments of the present invention, values for the characteristic parameters 
B and C in equation (8) are determined for blood bolus 22 using methods known in tiie art, 
such as for example a method described in Mourant et al; "Mechanisms of Light Scattering 
from Biological Cells Relevant to Noninvasive Optical-Tissue Diagnostics"; Applied Optics 
Vol 37. issue 16, pg 3586-3593, June 1998. In addition, a induced scattering coefficient 
a'sC^R). hereinafter a "reference scattering coefficient", for blood bolus 22 is determined for a 



12 



wo 2004/086965 



PCT/IL2004/000289 



reference wavelength A,r. In terms of the reference wavelength and associated reference 
scattering coefficient, the reduced scattering coefficient a'sCXj) may be expressed by 
(9) a's(Xi) = a's(XR)(XiaR)-C. 

In some embodiments of the present invention, a reference wavelength A,r for a tissue 
5 region is a wavelength for which the absorption coefficient aa(A,R) is known and the reference 
scattering coefficient a's(X-R) is determined from equation (2) and measurements of an 
extinction coefficient 0Ce(X,i0 at the reference wavelength. In some embodiments of the 
present invention, the absorption coefficient ocaCA^R) for a tissue region is known because the 
absorption coefficient of the tissue region is substantially determined by a component analyte 

10 whose concentration in the region is known. In some embodiments of the present invention, 
concentration of the component analyte is determined from a measurement of the extinction 
coefficient ae(Xi) for the region at a wavelength for which the extinction coefficient of the 
region is substantiaDy equal to the absorption coefficient of the component analyte. 
Optionally, as described above, measurements of the extinction coefficient ae(3lR) are 

15 acquired from time dependence of signals generated by transducers 42 responsive to 
photoacoustic waves stimulated in the region. 

By way of a numerical example, for detCTmining parameters of equation (9) required to 
determine a's(Xi) for blood, at 570 nm the magnitude of the reduced scattering coefficient 
a's(570) is between about 2 cm-1 and about Scm'l. The magnitude the absorption coefficient 

20 aa(570) of blood at 570 nm is about 280cm-l. The extinction coefficient oce(570) for blood at 
570 nm is therefore substantially equal to the absorption coefficient oca(570) of blood. In 
addition, the absorption coefficient of blood at 570 mn is substantially equal to tiie absorption 
coefficient of hemoglobin. Furthermore 570 nm is an isobestic wavelength for hemoglobin at 
which the absorption cross-sections for oxygenated and deoxygenated hemoglobin are about 

25 equal. Therefore, at 570 nm the concenti-ation of hemoglobin may be determined witiiout 
having to know the ratio of oxygenated hemoglobin to total hemoglobin from a measurement 
of the photoacoustic effect at 570 mn. Equation (4) for blood at wavelength 570 nm becomes, 
(10) oca(570) = ae(570) = Oah(570)xh, 

where Oah(570) is tiie absorption coefficient for hemoglobin at 570 nm and xh is flie 
30 concenti:ation of hemoglobin in blood. In accordance with an embodiment of tiie present 
equation (10) provides a value for ^q^. 
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810 nin is another isobestic wavelength in the absorption spectrum of hemoglobin at 
which the absorption coefficient of blood is also dominated by the absorption coefficient 
aah(810) of hemoglobin. However, at 810 nm the reduced scattering coefficient is not 
negligible and equation (2) becomes, 

(1 1) (Xe(810) = [3aah(810)xh(aah(810)xh + a's(810))] 1/2. 

Since is known from equation (10), equation (1 1) may be solved to provide a value, 
in accordance with an embodiment of the present invention, for the reduced scattering 
coefficient a's(810) at 810 nm. The scattering coefficient a's(^ at wavelength for blood 
may then be determmed by using 8 10 nm for the reference wavelength A,r and a's(810) for the 
reference scattering coefficient in equation (9) to provide, 

(12) a's(Xi) = a's(810) (V810)-C. 

To determine the coefficient C in equation (12), optionally, the extinction coefficient 
0Ce(A,) is determined for at least two other, non-isobestic, wavelengths of light at which 
hemoglobin concentration substantially determines the absorption coefficient of blood. 
Suitable wavelengths are preferably wavelengths that straddle 810 nm, for example 950 mn 
and 700 mn. Let the straddhng wavelengths be represented by X+ and A,", and collectively by 
X^. Let the ratio of oxygenated hemoglobin to total hemoglobin in the blood be represented by 
S and the absorption cross sections for oxygenated and deoxygenated at wavelengths be 
represented by OahoC^^ and OahDC^^ respectively, then the absorption coefficient OahHX^ 
for hemoglobm in the blood at wavelengths is, 

(13) cxah(X±)= [OahoC^'^S + aahD(A.^(l-S)]xh. 

Using equations (2), (12) and (13) and the extinction coefficients ae(X,i) determined 
for wavelengths A,± S and the exponent C may tbeaa. be determined from the two equations, 

(14) ae(A,±) = [3aah.(X±) (aah(X±) + a's(810)(^±/810)-C)]l/2. 

Substituting the right side of equation (9) for a'^iXO in equation (5) provides an 
equation of the form, 

(15) 2^o.iX.)x = l/2{a's(XR)(XiA,iO-C+ [a's(XR)2(Xi/A,R)-2C+(4/3)ae(Xi)2]l/2}. 

j 

Equation (15) for the Nx, mensuration wavelengths Xi provides a set of Ni 
simultaneous equations in the unknown concentrations xj, for each of which equations the 
right hand side the equation is known. In accordance with an embodiment of the invention, 
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controller 32 provides a value for the concentration xi of glucose responsive to constraints on 
the concentrations xi defined by the N^, simultaneous equations. Optionally, since water is a 
major component of living tissue and since the concentration of water is relatively labile at 
least one of the mensuration wavelengths is a wavelength, for example 1350nm, for which 
light is strongly absorbed by water and negligibly absorbed or scattered by other analytes in 
the body. 

For a number of mensuration wavelengths equal to a number Na of mensuration 
analytes, any of various well-known mettiods of manipulating and solving a set of 
simultaneous equations may be used to provide a value for xi. In some embodiments of the 
present invention a number of mensuration wavelengths N^, is greater than Na, resulting in a 
number of simultaneous equations greater than the Na unknown concentrations xj. For such 
cases a suitable best-fit algorithm, such as a least squares algorithm may be used to provide a 
solution for concentrations xj. 

In some embodiments of the present invention, equation (15) is treated as an equation 
in (Na+2) unknowns, where in addition to the unknown concentrations of the Na 
mensuration analytes, the reference scattering coefficient a's(A,R) and reference wavelength 
A,R are considered to be unknown constants. Measurements of the extinction coefficient 
0Ce(A^ are acquired for a plurality of N^ mensuration wavelengths Xi equal to or greater than 
(Na+2) to yield at least (Na+2) simultaneous equations of the form of equation (15). A set of 
at least (Na+2) simultaneous equation is sufficient to determine values for aU concentrations 
Xi, as weU as for a's(?lR) and A,r. Controller 32 provides a value for the concentration xi of 
glucose responsive to constraints on the concentrations xi, reference scattering coefficient 
«'s(^r) and reference wavelength Xr defined by the N;^ simultaneous equations. 

Similarly, ux accordance with some embodiments of tihe present invention, the 
exponent "C" in equation (15), is also considered to be an unknown and measurements of 
0Ce(Ai) are acqmred for at least (Na +3) mensuration wavelengths A^. The at least (Na+3) 
extinction coefficient measurements yield at least (Na+3) simultaneous equations of the form 
of equation (15), Controller 32 provides a value for the concentration xi of glucose responsive 
to constraints on the concentrations xf, reference scattering coefficient a's(A.R), reference 
wavelength A,r and exponent C defined by the N;^ = (Na+3) simultaneous equations. 

In some embodiments of the invention for which the scattering coefficient is expressed 
as an analytic fimction having at least one unknown characteristic parameter which is a 
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concentration of at least one of the mensuration analytes. OptionaUy, the concentration of at 
least one of the mensuration analytes includes tiie concentration of the target analyte. If the 
analytic function representing the scattering coefficient at wavelength X is written SCKPC), 
where X represents the set {xj} of concentrations of the mensuration analytes or a subset 
thereof, then equation (5) becomes, 

(16) 0Ca(Xi)= ^oJX.)x.= l/2{S(Ad,X)+[S(Xi,X)2+(4/3)ote(Xi)2]l/2}. 

j 

Similarly to the case of equation (5), the mensuration wavelengths provide a set of 
equations of the form of equation (16) in the Na unknown concentrations xj (l^j^ Na). 
The concentration of tiie target analyte is determined responsive to a solution of the set of 
equations. 

For some target analytes and conditions it is possible to choose an advantageous set of 
mensuration wavelength for determining concentration of an analyte in accordance with a set 
of equations of tiie form of equation (5) or equation 16. For example, as noted above at 570 
nm and 1350 mn the extinction coefficient for blood is substantially equal to the absorption 
coefficient of hemoglobin at 570 nm and water at 1350 nm respectively. At the isobestic 
wavelengtii 810 mn botii tiie absorption coefficient and tiie reduced scattering coefficient 
contribute to tiie extinction coefficient for blood. It is possible and can be advantageous to 
assay glucose, in accordance with an embodiment of flie invention, using tiiese tiiree 
wavelengths as mensuration wavelengths and hemoglobin, water and glucose as mensuration 
analytes. 

Iq particiUar, at 810 nm it can be advantageous to use a set of simultaneous equation at 
tiie mensuration wavelengtiis for which at least one equation has tiie form of equation (16) and 
tiie reduced scattering coefficient is represented by an analytic function S(Xi,X). For example, 
if tiie concentrations of hemoglobin, water and glucose are represented by xh, x^ and Xg, 
respectively, S(^,X) in equation (16) optionally becomes S(^,xh,XwXg). OptionaUy, 
S(^,xh,x^ Xg) may be expanded in a Taylor series to a desired order in the concentrations jqj, 
x^ and Xg, Coefficients in the Taylor series may be determined from a suitable model and/or 
empirically. For example, the coefficients may be determined using an expression for the 
reduced scattering coefficient described in "Dynamic optical coherence tomography in studies 
of optical clearing, sedimentation, and aggregation of immersed blood"; Valery V. Tuchin, 
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Xiangqun Xu, and Ruikang K. Wang; APPLffiD OPTICS Vol. 41, No. 1 , 258-271, January 
2002. Optionally, since at wavelengJhs 570 nm and 1350 mn the extinction coefficient is 
dominated by the absorption coefficient the reduced scattering coefficient is assumed to be 
zero and an expression for S(^,xii,Xw Xg) is used only in an equation of the form (16) at 810 



nm. 



In the above examples, extinction coefficients ae(Ai) for the mensuration 
wavelengths and for the reference wavelength that are used to determine glucose 
concentration xi are described as being determined using the photoacoustic effect. In some 
embodunents of the present invention optical coherence tomography (OCT) is used to 
determine at least one of the extinction coefficients used to determine concentration of an 
analyte. 

OCT generally provides signals for determining an extinction coefficient having better 
SNR than photoacoustic effect signals at wavelengths for which the extinction coefficient is 
determined substantiaUy by a reduced scattering coefficient Photoacoustic effect signals 
generaUy have better SNR than OCT signals at wavelengths for which an extinction 
coefficient is dominated by an absorption coefficient In accordance with an embodiment of 
the present invention, a glucometer for assaying glucose in a tissue region comprises at least 
one acoustic tiransducer and m addition an "OCT" interferometer. Photoacoustic signals 
generated by the at least one acoustic transducer are processed to determine extinction 
coefficients used to assay glucose for wavelengths at which an absorption coefficient 
dominates in determining a value for the extinction coefficient Interference signals generated 
by the interferometer are processed to determine extinction coefficients used to assay glucose 
for wavelengths at which a scattering coefficient dominates in determining a value for tiie 
extinction coefficient 

Fig. 2 schematicaUy shows a glucometer 100, in accordance with an embodiment of 
the present invention comprising at least one acoustic transducer and an OCT interferometer. 
The components of the OCT interferometer are shown in a very schematic and simpHfied 
manner. Glucometer 100 is schematically shown assaying glucose in blood bolus 22 in blood 
vessel 23 of tissue region 26. 

Glucometer 100 optionaUy comprises a controller 102, and a hght source 104, 
optionally located in tiie controller, that provides semi-coherent Ught at wavelengths for which 
it is desired to determine an extinction coefficient to assay glucose, in accordance with tiie 
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invention. An optic fiber 36 is coupled to Ught source 104 via an optical coupler 106. An end 
38 of fiber 36 is optionally mounted to a support structure 40 to which acoustic transducers 42 
are mounted. 

To assay glucose in bolus 22 controller 102 controls light source 104 to transmit at 
least one pulse of Ught into optical fiber 36 at each of wavelength for which it is desired to 
determine an extinction coefficient for bolus 22. An optical coupler 108 couples a portion of 
Ught transmitted by Ught source 104 along optic fiber 36 to an optical fiber 1 10 and transmits 
a portion of the light towards end 38 of fiber 36 firom which end the light exits the fiber to 
illuminate tissue region 24. 

Some of the Ught that is transmitted along fiber 36 to exit the fiber at end 38 is 
absorbed in tissue region 26 and stimulates photoacoustic waves 52 in the region and some of 
the light is scattered by material in region 26. As in glucometer 20 acoustic transducers 42 
generate signals responsive to photoacoustic waves 52. Some of the scattered Ught reenters 
fiber 36 at end 38 and propagates back towards controUer 102 through optical coupler 108. 

Light coupled to optical fiber 1 10 by coupler 108. exits the fiber from an end 1 12 and 
is reflected back into the fiber by a mirror 1 14. A portion of the light reflected back into optic 
fiber 1 10 is directed by optical coupler 108 to controller 102. When Ught reflected from miiror 
1 14 and scattered Ught from tissue region 26 that reenters optic fiber 36 reaches controller 
102, the light is directed by coupler 106 to a combiner 116. Combiner 116 superposes the 
scattered light and the reflected Ught at an interference region (not shown) to generate an 
interference signal. The position of mirror 114 relative to fiber end 112 is controlled by 
controUer 102 to determine a desired path length fix)m Ught source 104 to mirror 1 14 and back 
to combiner 1 16. The path length is determined so that substantially only light that is scattered 
in tissue region 26 from desired locations in the region generates an interference signal. Semi- 
coherent Ught source 104, couplers 106 and 108, optic fibers 36 and 110, mirror 114 and 
combine 1 16 cooperate to function as an OCT interferometer. 

ControUer 102 processes signals generated by acoustic transducers 42 to determine 
which of the signals correspond to photoacoustic waves originating at different locations in 
bolus 22. ControUer 102 controls the position of mirror 114 to scan bolus 22 and generate 
interference signals corresponding to Ught reflected fix>m different locations in the bolus. 

For Ught at mensuration wavelengths having an extinction coefficient dominated by an 
absorption coefficient controUer 102 optionaUy uses signals generated by transducers 42 
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corresponding to locations in bolus 22 to detennine an extinction coefficient at the wavelength 
for the bolus. For light at mensuration wavelengths having an extinction coefficient dominated 
by a scattering coefficient, controller 102 optionally uses interference signals generated by 
combiner 116 corresponding to locations in bolus 22 to detennine an extinction coefficient at 
the wavelength for the bolus. The controller uses extinction coefficients to assay glucose in 
the same way that glucometer 20 uses extinction coefficients to assay glucose. 

In the description and claims of the application, each of the verbs, "comprise" 
"include" and "have", and conjugates thereof, are used to indicate that the object or objects of 
the verb are not necessarily a complete listing of members, components, elements or parts of 
the subject or subjects of the verb. 

The present invention has been described using detailed descriptions of embodiments 
thereof that are provided by way of example and are not intended to limit the scope of the 
invention. The described embodiments comprise different features, not all of which are 
required in all embodiments of the invention. Some embodiments of the present invention 
utilize only some of the features or possible combinations of the features. Variations of 
embodiments of the present invention that are described and embodiments of the present 
invention comprising different combinations of features noted in the described embodiments 
will occur to persons of the art. The scope of the invention is limited only by the following 
claims. 
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